Dynamics of nearly spherical vesicles in an external flow 
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We analytically derive an equation describing vesicle evolution in a fluid where some stationary flow is ex- 
cited regarding that the vesicle shape is close to a sphere. A character of the evolution is governed by two 
dimensionless parameters, 5* and A, depending on the vesicle excess area, viscosity contrast, membrane vis- 
cosity, strength of the flow, bending module, and ratio of the elongation and rotation components of the flow. 
We establish the "phase diagram" of the system on the S — A plane: we find curves corresponding to the tank- 
treading to tumbling transition (described by the saddle-node bifurcation) and to the tank-treading to trembling 
transition (described by the Hopf bifurcation). 

PACS numbers: 87.16.Dg, 87.17.Jj, 05.45.a 



Vesicles are closed membranes which separate two regions 
occupied by possibly different fluids. The vesicles are attract- 
ing significant attention not only due to their resemblance with 
biological objects but also because of their importance in dif- 
ferent industries such as pharmaceutics where they are used 
for drug transportation. A natural problem which arises in 
these applications is understanding of how a single vesicle 
behaves in an external flow. This non-equilibrium problem 
has revealed a variety of new physical effects and became a 
subject of intense experimental and theoretical studies. Lab- 
oratory experiments Uj, |2l |3l |4|, jSJ have shown that the vesi- 
cles immersed in a shear flow exhibit at least two qualitatively 
different types of behavior, either tank-treading or tumbling 
motion. In the tank-treading regime a vesicle shape is station- 
ary, it is ellipsoid oriented at an angle with respect to the shear 
flow. In the tumbling regime the vesicle experiences periodic 
flipping in the shear plane. A novel type of behavior: trem- 
bling, discovered in the work |0] is an intermediate regime 
between tank-treading and tumbling in which a vesicle trem- 
bles around the flow direction. 

Constructing a phase diagram for all these regimes depend- 
ing on the external parameters is a challenging and an ex- 
tremely difficult task because the problem in consideration 
is both strongly non-linear and non-equilibrium. As long as 
no analytic solution of this problem exists theoretical stud- 
ies were based either on numerical simulations or on some 
approximations allowing analytical treatment. Numerical in- 
vestigations of this problem involved several different com- 
putational schemes, including boundary element method j?!, 
mesoscopic particle-based approximation |^ [13, 11. 121. 
and an advected field approach 11311 . These approaches have 
shown qualitative agreement with experiments however did 
not solve the problem of constructing the vesicle dynamics 
phase diagram completely. Analytical studies of the prob lem 
can be divided in two major classes. In the first one ll9[ ll4l[l5|] . 
phenomenological models of a vesicle dynamics based on the 
classical work of Keller and Skallak lll6ll were proposed and 
proved themselves to be rather efficient in explaining the ex- 
periments. In the second series of works ilTj 18, 19, 20] the 
studies focused on quasi-spherical vesicles whose shape can 
be parameterized by spherical harmonics expansion. A pertur- 



bation scheme around the Lamb solution for spherical particle 
in external flow allows one to derive the dynamic equations 
for the shape and the orientation of a vesicle and to inves- 
tigate them analytically. In this letter we propose a natural 
extension of the theory developed for quasi-spherical vesicles 
by accounting higher-order expansion terms. We show that 
these additions produce a qualitative change in the phase di- 
agram and make it significantly more complicated. The re- 
sulting diagram which is the main result of this paper contains 
all three types of vesicle behavior which were observed ex- 
perimentally. In our work we analyze how the vesicle dynam- 
ics depends on different control parameters, such as viscos- 
ity contrast, vesicle excess area, internal membrane viscosity, 
strength of the flow, and ratio of the elongation and rotation 
components of the flow. We analyze also the vesicle orienta- 
tion in the tank treading regime. 

Speaking about membranes we have in mind lipid bilay- 
ers. Physical properties of such objects have been extensively 
studied both experimentally and theoretically (see e.g. the 



books 112 IL 1221 12311 and the reviews 112411 ). There are several 



features of the membranes which are important for our anal- 
ysis. First, we assume that the membrane is in a fluid state 
(is a 2d liquid), which is typical of the lipid bilayers under 
normal conditions. Second, we assume that the vesicle has an 
excess area, that enables one to treat the vesicle membrane as 
incompressible. Third, we assume that the membrane is im- 
permeable to the surrounding liquids, the condition is usually 
well satisfied in experiment. Finally, we take into account the 
membrane internal viscosity, which could play an essential 
role, say, in a vicinity of the lipid-bilayer melting point ll25ll 
(see alsoRefs. Hill). 

We assume that the vesicle size is much larger than the 
membrane thickness. Then, in the main approximation, the 
vesicle membrane can be treated as infinitesimally thin, that is 
as a 2c? object immersed into a 3d fluid. By other words, in the 
limit the vesicle can be considered as an interface separating 
two generally different fluids. The two properties, the mem- 
brane incompressibility and impermeability, imply that both 
the vesicle volume V and the membrane area A are conserved. 
The excess area can be characterized by a dimensionless fac- 
tor A, which is defined as ^ = {Att + A)rQ, where ro is a 
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vesicle "radius" determined by its volume: V = 47rrQ/3. The 
excess area is non-negative, A > 0, and the minimal value 
A = corresponds to ideal spherical geometry. 

The membrane free energy can be written as the following 
integral over the membrane position 



dA[G + ^H^ 



(1) 



where a is the membrane surface tension, H is its mean cur- 
vature, and K is Helfrich module. The last term in Eq. ([T]i 
describes energy related to membrane bending distortions. 
Note that the surface tension <t is a quantity adjusting to other 
membrane parameters (similarly to pressure in incompressible 
fluid) to ensure a given value of the membrane area A. 

The membrane moves together with surrounding fluid that 
is the velocity field v is continuous on the membrane and v 
determines the membrane velocity as well as the fluid veloc- 
ity. We divide the flow near the vesicle into two parts: an 
external flow which would occur in the fluid in the absence 
of the vesicle, and an induced flow which is excited as a re- 
sult of the vesicle reaction to the external flow. We assume 
that a characteristic scale of the external flow is much larger 
than the vesicle size. Then the external flow velocity V near 
the vesicle can be approximated by a linear profile. Generally, 
the external flow has both the strain (elongation) and the ro- 
tational contributions: diVk = Sik + eikji^j, where s is the 
(symmetric) strain matrix and w is the angular velocity vec- 
tor. The strain can be characterized by its strength s, defined 
as = tr Note that for a shear flow s ^ lo = 7/2, 

where 7 is the shear rate. 

We examine nearly spherical vesicles, that is the excess 
area parameter A is regarded to be small. Then it is nat- 
ural to describe the vesicle shape (membrane position) as 
r = ro[l + <p)] where r, 0, </3 are spherical coordinates in 
the reference system with origin chosen at the vesicle center. 
The quantity u is a dimensionless displacement characteriz- 
ing deviations of the membrane shape from a spherical one. 
For small A the value of u can be estimated as w \/A. 
Therefore in our scheme u ^ \ and one can formulate a per- 
turbation theory in the parameter. 

We assume that both interior fluid and exterior one are 
Newtonian and that Reynolds number is vanishingly small, so 
the fluids can be described in terms of Stokes equation. We as- 
sume also that an adiabaticity condition max{s, k/ (r/rp)} <C 
77/ (prg) is satisfied where ry is viscosity and g is mass density 
of the exterior fluid. The same condition is assumed for the 
interior fluid. Then the flow inside and outside the vesicle can 
be treated as instantaneously adjusted to the vesicle motion. 
In this case it is possible to establish a closed dynamic equa- 
tion for the membrane displacement u. For the purpose one 
should find the velocity field v inside and outside the vesi- 
cle at a given displacement u{9, (p) and then equate dfU to 
the membrane normal velocity. To find the velocity field one 
should solve the Stokes equation with boundary conditions 
dictated by the membrane incompressibility and by the mo- 
mentum balance that includes membrane forces determined 



by the energy 

(IB EESE 1|] and membrane viscosity. To 
realize the program for a nearly spherical vesicle one can use 
a generalization of the Lamb scheme applicable to a spheri- 
cal body, see Ref. ^2^- As a result, one finds the dynamic 
equation for u as a series in u. 

In the main approximation in u one obtains 



1 <5.F(3) 



Su 
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where a is some dimensionless operator, reflecting all viscous 
mechanisms, the Z-axis of our reference frame is chosen to be 
directed opposite to the angular velocity u), and JF^-^' is an ex- 
pansion up to third order in u of the free energy ([T])- Note that 
the elongation and the rotation parts of the external flow are 
separated in Eq. the angular velocity u; extends the time 
derivative whereas the strain determines the non-uniform term 
playing a role parallel to the free energy derivative. For an ex- 
ternal shear flow SijriVj/r^ = (7/2)sin^6' sin(2(/3), where 
X-axis is directed along the velocity. 

The two differences between the equation ^ and the analo- 
gous equations obtained by Misbah 11911 and Vlahovska et. al. 
12011 are the inclusion of membrane viscosity and third order 
expansion term of Helfrich energy. An account of the third or- 
der term changes the phase diagram completely since the sec- 
ond order term becomes vanishingly small (due to the surface 
tension adjustment) in the vicinity of the tumbling to tank- 
treading transition region predicted in Refs. 20D- An im- 
portance of the third-order term is discussed also by Noguchi 
and Gompper 115 1 . 

We take into account only the main contribution to u deter- 
mined by second order angular harmonics. It is justified by 
the smallness of A and by the fact that the non-uniform term 
in the right-hand side of Eq. (|2]i is a linear combination of 
second order angular harmonics. The operator a in the case is 
reduced to a constant 
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where fj is viscosity of the interior fluid and C, is membrane 
viscosity. Note that the expression ([3]) includes the viscosity 
contrast 7^/77. After passing to the variable u/^fK the equa- 
tion (|2]i acquires a self-similar form containing the parame- 
ters ^/Aauj/s, s7]rQ/{KA) and some dimensionless param- 
eters characterizing ratios of the eigenvalues of the matrix s 
and mutual orientation of the vorticity vector uj and the main 
axes of the matrix s. 

Below, we consider a plain flow where the external velocity 
V lies inX — Y plane. In the case the equation for u/\/A de- 
pends on two dimensionless parameters which can be chosen 



A = 



x/3 
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In the case of weak strains, S <^ 1, the vesicle conserves its 
equilibrium shape, and in the case of strong strains, S ^ 1, 
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FIG. 1: Phase diagram on the S — A plane. 



the vesicle shape is determined by the external flow. The pa- 
rameter A determines an effectiveness of the rotational part of 
the external flow. Note that our theory is applicable, partic- 
ularly, to purely elongation flow, where u; — and, conse- 
quently, A = 0. 

For the plain flow, one can further reduce description of the 
vesicle shape which is dependent on two parameters, and 
$: 



V5A 
4^/2^ 



sinO 



(1-3 cos^e) +cose sin^6' cos2((^- $) 



Here, the axes X and Y are chosen to be directed along the 
principal axes of the strain matrix s. Note that the "angles" 
$, 6 determine the main axis direction of the vesicle projec- 
tion to the X—Y axis. In the limit A <C 1 the inclination angle 
of the vesicle projection main axis is (ys* = $ for cos > 
and if* — $ + 7r/2 otherwise. The parameter Q determines 
also the vesicle shape and can be directly related to the pa- 
rameter R = cos 8 introduced in Ref. 11911 . In terms of the 
variables $ and 0, Eq. (|2]l is reduced to a couple of equations 



rdt^ = ^ 



cos(2$) 
cos 9 



A 



TdtQ = -S'sinesin(2$) +cos(3e), 
wncrc T — 
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The tank treading vesicle motion corresponds to stable sta- 
tionary points of the equations ( |5l6t . Equating to zero the right 
hand sides of the equations one finds relations determining a 
stationary point at given parameters S and A. To investigate 
stability of the point one should linearize the equations (|5] 
|6ll neai- the point to obtain rdtiSQ, 6^) = B{SO, (5$). The 
point is stable, if both eigenvalues of the matrix B have nega- 
tive real parts. Thus the stability conditions are tr S < and 
dot B > 0. Regions in the S* — A plain where stable points 



exist are indicated in Fig. [Dby blue color. Note that there ex- 
ists an additional region of stationary points which are stable 
in terms of and <!'. However, a stability investigation in the 
framework of Eq. (|2|i shows that the points are unstable in the 
extended space where the spherical harmonics 1^2, ±i 0) are 
involved into play. Therefore we do not analyze the points. 

If A is increased at a given value of S then first the sta- 
ble point is in the domain where $ > 0. At the boundary of 
the domain, at A = 2/\/3, the parameters $ and come 
to the values <I> = 0, = tt/6 . Then the stable point 
moves into the domain, corresponding to A > 2/V3 and 
$ < 0. The green line in Fig. [T]is determined by the condition 
{tr B)^ = 4deti3, it separates regions where oscillatory and 
damping regimes of approaching the stable point occur The 
boundary of the stability domain can be determined either by 
the condition det-B = or tri? = 0. The corresponding 
lines are depicted in Fig. [Uby yellow and red lines, respec- 
tively. The first case is realized at small S whereas the second 
case is realized at large S. At increasing A a Hopf bifurca- 
tion occurs at the red line whereas a saddle-node bifurcation 
occurs at the yellow line. A limit cycle is realized in the sys- 
tem for points above the red and the yellow lines (that is at 
large A). Following the work |0] we distinguish tumbling and 
trembling regimes. In our terms, trembling corresponds to 
a limit cycle where (p varies in a restricted interval whereas 
tumbling means an increase of ip by 2tt during the cycle. The 
dashed blue line in Fig. [T] separates the two regimes, trem- 
bling is realized below the line. A possibility of such regime 
is discussed also by Misbah I19I1 . The sequence of the vesicle 
states predicted by our theory corresponds to one observed in 
experiment Qualitative the picture is also similar to the 
one observed in recent numerical simulations by Noguchi and 
Gompper lIlSll . Note that at infinitely large fj ot ( a solid ball 
behavior of the vesicle should be observed which is tumbling 
as Jeffery established |33]). In our scheme large ^ or C implies 
large A, corresponding to tumbling region, indeed. 

A special point S = \/3, A = 2/\/3, where two transi- 
tion lines terminate, is analogous to some extend to fluid crit- 
ical point. Say, frequency of the Hopf bifurcation behaves 
cx \/<5'^ ~ 3. Similar "smoothing" of the saddle-node bifurca- 
tion near the point occurs. A complicated phase structure is 
observed near the special point, its narrow vicinity is plotted 
in the insertion to Fig. [T] The green region corresponds to 
coexistence of two stable points. There exists also a region 
where the stable point and the limit cycle coexist. The limit 
cycle becomes unstable at the dashed green line. 

Thermal fluctuations of the vesicle shape can be investi- 
gated in spirit of Refs. IhI HS (see also the book Q). 
The fluctuations produce isotropic Brownian motion on the 
sphere, defined by the "angles" and $ with the character- 
istic angular diffusion coefficient D ^ kT/ (jyrQ A), where T 
is temperature. In the tank-treading regime and strong exter- 
nal flow, 5^1, fluctuations induced by this diffusion can 
be estimated as {5Qf - {5^f - Dt/S if A < 2/\/3 
and as {SQf - {5^f ~ Dr if A > 2/^3, where r is 
defined by Eq. (|7]i. For weak external flow, 5* <C 1, fluctu- 
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ations of vesicle shape are given by estimate (SQ)"^ ^ Dt, 
whereas vesicle orientation experiences more stronger fluc- 
tuations, (S^)^ ^ Dt/S. The thermal fluctuations play an 
essential role near the tank-treading to tumbling (or tank- 
treading to trembling) transition "smearing" it. Results of 
detailed investigation of this phenomenon will be published 
elsewhere. 

To conclude, we have theoretically investigated the vesicle 
dynamics in an external stationary flow. The general scheme 
based on solving the 'id hydrodynamic (Stokes) equations 
with boundary conditions posed on the membrane enabled 
us to analyze in detail dynamical properties of nearly spher- 
ical vesicles. We constructed the phase diagram of the sys- 
tem in terms of two dimensionless parameters which are 
combinations of physically observed quantities. We demon- 
strated that there are two different regimes realized in weak 
and strong external flows. In weak flows the vesicle shape is 
close to an equilibrium one which is a prolate ellipsoid and a 
role of the external flow is reduced mainly to an orientation of 
the ellipsoid. Then the tank-treading to tumbling transition oc- 
curs which can be described by a saddle-node bifurcation. In 
strong flows the vesicle shape and orientation are determined 
by the flow. Then the tank-treading to trembling transition oc- 
curs which can be described by a Hopf bifurcation. The two 
transition curves are separated by a point vicinity of which 
needs a special consideration. 
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